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ABSTRACT 

We investigate the kinematic signatures induced by spiral and bar structure in a set of 
simulations of Milky Way-sized spiral disc galaxies. The set includes test particle sim¬ 
ulations that follow a quasi-stationary density wave-like scenario with rigidly rotating 
spiral arms, and A^-body simulations that host a bar and transient, co-rotating spiral 
arms. From a location similar to that of the Sun, we calculate the radial, tangential 
and line-of-sight peculiar velocity fields of a patch of the disc and quantify the fluctu¬ 
ations by computing the power spectrum from a two-dimensional Fourier transform. 
We find that the peculiar velocity power spectrum of the simulation with a bar and 
transient, co-rotating spiral arms fits very well to that of APOGEE red clump star 
data, while the quasi-stationary density wave spiral model without a bar does not. We 
determine that the power spectrum is sensitive to the number of spiral arms, spiral 
arm pitch angle and position with respect to the spiral arm. However, it is necessary 
to go beyond the line of sight velocity field in order to distinguish fully between the 
various spiral models with this method. We compute the power spectrum for different 
regions of the spiral discs, and discuss the application of this analysis technique to 
external galaxies. 

Key words: galaxies: evolution - galaxies: kinematics and dynamics - galaxies: spiral 
- galaxies: structure 


1 INTRODUCTION 

Understanding how large disc galaxies such as the Milky 
Way evolved into their present state is one of the principal 
aims of galactic astrophysics. In the Milky Way, particu¬ 
larly well known constraints for galaxy evolution include, 
for example, the correlations of the age of stars with their 
kinematic and chemical properties, i.e., the age-velocity dis¬ 
persion and age-metallicity relations. Holmberg et al. (2009) 
showed from a combination of radial velocities derived from 
the metallicity data of F- and G-dwarf stars (Nordstrom 
et al. 2004) and distances and proper motions from Hippar- 
cos (van Leeuwen et al. 2007), that the velocity dispersion 
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of stars increases with stellar age in all directions. Although 
the determinations of stellar age can be uncertain (e.g. Reid 
et al. 2007; Casagrande et al. 2011), this positive trend be¬ 
tween age and velocity dispersion is supported by other stud¬ 
ies (e.g. Dehnen & Binney 1998; Seabroke & Gilmore 2007; 
Aumer & Binney 2009). Among others, a possible explana¬ 
tion for these trends is that a combination of spiral structure 
(Barbanis & Woltjer 1967; Carlberg & Sellwood 1985) and 
giant molecular clouds (Spitzer & Schwarzschild 1951) ef¬ 
fectively scatter stars in the planar and vertical directions 
respectively, thus increasing the velocity dispersion of the 
stars while maintaining the observed radial to vertical ve¬ 
locity dispersion ratio. However, the mechanisms that drive 
these trends are still unknown. 

Although the observed stellar kinematics are often de¬ 
scribed by a symmetric Gaussian in the radial and vertical 
directions, and a skewed Gaussian in the rotational direc- 
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tion, there is evidence that the velocity distribution in planar 
velocity space exhibits multiple streaming structures (e.g. 
Dehnen & Binney 1998; Chereul et al. 1998, 1999; Famaey 
et al. 2005; Antoja et al. 2008; Bovy et al. 2009; Antoja 
et al. 2015), which have been attributed to the gravitational 
perturbations provided by non-axisymmetric structure. In 
particular, it has been reported that the local moving group 
known as the ‘Hercules’ stream is likely caused by the con¬ 
stant periodic perturbations supplied by the outer Lindblad 
resonance of the bar (e.g., Dehnen 2000; Fux 2001; Minchev 
et al. 2010), which is thought to be close to the solar radius. 
The connection between the bar and such a moving group 
has been used to pin down some bar parameters, such as 
the pattern speed and orientation angle with respect to the 
solar position, either from the determination of the moving 
group position (e.g., Dehnen 2000; Bovy 2010; Monari et al. 
2013; Antoja et al. 2014) or the variation of the local Oort 
constants (Minchev & Quillen 2007). 

Aside from the bar, spiral arms have been highlighted 
also to be the source of non-circular streaming motions (see 
for example De Simone et al. 2004; Bovy & Hogg 2010), 
and in some cases has been directly linked to local moving 
groups such as the ‘Hyades’ moving group (e.g., Sellwood 
2010; Pompeia et al. 2011). In particular, spiral structure 
has been shown to change the guiding radii of disc stars by a 
process known as radial migration (Lynden-Bell & Kalnajs 
1972; Sellwood & Binney 2002). This occurs near the co¬ 
rotation radius of spiral arms, where the tangential gravita¬ 
tional force from spiral arms introduces a steady torque that 
imparts angular momentum changes on stars. Moreover, ra¬ 
dial migration may be amplified if co-rotation resonances of 
spiral arms overlap other resonances such as the bar outer 
Lindblad resonance (e.g., Minchev & Famaey 2010). Because 
radial migration occurs around regions of co-rotation, its ef¬ 
fects depend on the location and number of resonances, and 
therefore the nature of the spiral arm itself. 

The most widely accepted theory of spiral structure 
is quasi-stationary density wave theory (Lin & Shu 1964; 
Berlin et al. 1989), in which the spiral arm is described as a 
stationary wave that rotates around the galaxy with a con¬ 
stant pattern speed without changing shape. However, there 
are conflicting results in the literature regarding the valid¬ 
ity of density wave theory in external galaxies. For example, 
application of the Tremaine-Weinberg method (Tremaine & 
Weinberg 1984) to external inclined galaxy discs typically 
favours radially decreasing spiral pattern speeds (e.g., Mer- 
rifield et al. 2006; Speights & Westpfahl 2011). Furthermore, 
Foyle et al. (2011) and Ferreras et al. (2012) analysed the 
distribution of star forming tracers across spiral arms in ex¬ 
ternal galaxies and found no evidence of the angular offsets 
between different star forming tracers (although see Egusa 
et al. 2009) predicted by theories of long-lived, rigidly ro¬ 
tating spiral arms (Roberts 1969). From a numerical per¬ 
spective, Dobbs & Bonnell (2008) and Dobbs et al. (2010) 
performed numerical simulations of isolated spiral galaxies 
and confirmed that tracer offsets are present in simulations 
with an imposed rigidly rotating spiral arm potential, like 
density waves, while they are absent in those without (see 
also Grand et al. 2012b). However, some numerical studies 
of cosmologically simulated disc galaxies exhibit a distribu¬ 
tion of young stars consistent with the predictions of classic 
density wave theory for long-lived spirals (Pilkington et al. 


2013), and the velocity gradient found in the RAVE sur¬ 
vey (Siebert et al. 2011) has been shown to be reproducible 
with Lin-Shu density waves (Siebert et al. 2012). On the 
other hand, it is not clear that this is also the case for ve¬ 
locity fluctuations on larger scales, which we will consider in 
the present paper. 

In general, spiral arms in numerical simulations do not 
reproduce the long-lived single density wave structure (e.g., 
Sellwood 2011), and their scrutiny has spurred alternative 
ideas of the nature of spiral structure, such as multiple 
mode theory (Minchev et al. 2012; Comparetta & Quillen 
2012; Sellwood & Carlberg 2014), manifold theory (Romero- 
Gomez et al. 2007; Athanassoula et al. 2009) and non¬ 
steady, co-rotating spiral arms (Wada et al. 2011; Grand 
et al. 2012a,b; Baba et al. 2013). In the last description, 
spiral arms are transient, recurrent features that wind up 
and disappear with time, continually being replaced by new 
transient features that form from small perturbations ampli¬ 
fied by a mechanism akin to swing amplification (Goldreich 
& Lynden-Bell 1965; Julian & Toomre 1966; Toomre 1981; 
D’Onghia et al. 2013). These structures are found to be ap¬ 
proximately co-rotating with the disc stars at every radius 
(e.g.. Grand et al. 2012a,b), and cause continual radial mi¬ 
gration at every radius that manifests as large systematic 
streaming motions along the spiral arm (Grand et al. 2015; 
Hunt et al. 2015), which is expected to be different from the 
peculiar motions induced by density wave spirals. The pecu¬ 
liar motions associated with different spiral arm models may 
therefore offer distinguishing observational predictions that 
can be tested particularly well in the Milky Way - the only 
galaxy for which we are able to get star by star informa¬ 
tion. The effects of different spiral arm models on the veloc¬ 
ity distribution must be combined with those of spiral arm 
parameters such as the number of spiral arms, constraints 
on which have already been placed by analysis of peculiar 
motions of gas (e.g. Pettitt et al. 2014) and stars in pencil 
beam surveys (e.g. Minchev & Quillen 2008) in numerical 
simulations. 

A useful way to characterise the peculiar motions in 
the disc was introduced by Bovy et al. (2015), who stud¬ 
ied the kinematics of red clump stars from APOGEE (Al- 
lende Prieto et al. 2008), abbreviated APOGEE-RG data, 
and computed the power spectrum of the peculiar line-of- 
sight (LOS) velocity field. From a comparison between the 
observed power spectrum and that of various disk galaxy 
models, they concluded that the data favoured a bar only 
model that produced a peak peculiar velocity of 12 km s“^ 
on scales of 2 kpc. However, a successful model for the 
Milky Way ought to take into account the bar and spiral 
structure, both of which imprint their respective kinemat- 
ical signatures in the peculiar velocity field. It is therefore 
key to model the perturbative effects of different spiral arm 
models both with and without a bar in order to understand 
the individual effects of each component. In this paper, we 
aim to build on the power spectrum analysis performed by 
Bovy et al. (2015) in order to compare the peculiar veloc¬ 
ity fields produced in a range of spiral arm models with the 
APOGEE-RC LOS velocity data. In particular, we seek to 
establish whether or not spiral arms that form in A^-body 
simulations are able to reproduce the characteristic features 
of the peculiar LOS velocity power spectrum. In addition, 
we decompose the peculiar velocity field into radial and tan- 
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gential components. Scrutiny of these velocity fields allows 
us to quantify the induced peculiar motions in each model, 
in particular the A'^-body and spiral density wave models, 
and make observational predictions of each one that may 
be tested with Gaia (Perryman et al. 2001; Lindegren et al. 
2008) in the near future. Furthermore, the simulations al¬ 
low us to probe the effect of spiral arm pitch angle, spiral 
arm number and solar position with respect to the spiral 
arm on the resulting power spectrum, which may provide 
constraints on spiral arm parameters. We take further ad¬ 
vantage of the simulations by computing the power spectra 
of both stars and gas over the whole disc (not just the solar 
neighbourhood region), which serve as observational predic¬ 
tions of the spiral-bar induced kinematics for external galax¬ 
ies that may be tested by integral field spectroscopy (IFU) 
surveys such as MUSE (Bacon et al. 2010) and CALIFA 
(Sanchez et al. 2012). 

We note here that a limitation of this work is that we 
do not consider any external perturbations such as merg¬ 
ers and satellite interactions. However, satellite-disc cross¬ 
ings have been shown to generate, at a given position in 
the galaxy disc, ring-like features in planar velocity space 
(Minchev et al. 2009; Gomez et al. 2012) and waves in ver¬ 
tical velocity (Widrow et al. 2012; Gomez et al. 2013). It 
is therefore possible that such interactions can contribute 
to the observed peculiar velocity features that we discuss in 
this paper. However, an investigation into this possibility is 
beyond the scope of this paper, and we leave this to future 
studies. 


2 SIMULATIONS 

In this paper we analyse five simulations of isolated Milky 
Way-sized disc galaxies; three A-body simulations and two 
test particle simulations. The A-body simulations comprise 
a live disc component embedded in a static dark matter halo 
potential that follows the NFW profile (Navarro et al. 1997). 
The dark matter halo is modelled as a static potential in 
order to reduce computational costs and to avoid scattering 
of disc particles by dark matter particles, which are often 
restricted by computational limitations to be much more 
massive than disc particles. Further details of the simulation 
code are available in Kawata & Gibson (2003); Kawata et al. 
(2013, 2014). 

The main simulation parameters relevant for this study 
are summarised in Table 1. The fiducial simulation, labelled 
K14 (see Kawata et al. 2014, for more details), develops a 
bar and a prominent two-armed spiral structure. The spi¬ 
rals are transient features that wind up and disrupt owing 
to their co-rotation with the disc material. However, they 
are also recurrent features, and therefore the two-armed spi¬ 
ral morphology is sustained over the simulation duration. 
Simulations G13F (simulation F presented in Grand et al. 
2013) and G14 (presented in Grand et al. 2014) develop 
transient spiral structure only, i.e., there is no bar. A fur¬ 
ther difference between these simulations is that simulation 
G13F has a central static stellar bulge component of mass 
Mb,-k = 4.0 X 10^° Mq, and the dark matter distribution is 
slightly less centrally concentrated than in simulation G14. 
The difference in halo concentration means a lower disc-halo 
mass fraction at radii A < 10 kpc in simulation G14 in com¬ 
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Figure 1. Two-dimensional schematic of the coordinates systems 
and peculiar velocity definitions described in the text. The Galac¬ 
tic centre and Solar position are marked by the solid circle and 
encircled dot, respectively. The line-of-sight to a star at solar ra¬ 
dius in the disc plane and galactic longitude, I, is represented by 
the dashed line. The direction of rotation is clockwise. 


parison to G13F, which leads to a higher number of spiral 
arms in the former. 

To determine the spiral arm strength in each model, we 
calculate the radial force per unit mass, F_r, at a series of 
points distributed evenly in azimuth in a ring at a Galacto- 
centric radius of 8 kpc, and then the ratio 


q{R — 8 kpc) 


\Fn{e)-Fn\me.^ 


Fr 


R=8 kpc 


( 1 ) 


which describes the amplitude of radial force fluctuation 
caused by non-axisymmetric structure, |Fr( 6') -F itlmax, as 
a fraction of the total axisymmetric radial force, Fr. We 
refer to q as the radial force ratio, and take it to be the indi¬ 
cator of the strength of non-axisymmetric structure in each 
simulation. 

Simulation F14, presented in Fame et al. (2014), is a 
test particle simulation in which the orbits of massless disc 
particles are integrated in a superposition of background 
axisymmetric (disc-|-bulge-|-dark matter halo) and spiral po¬ 
tentials. The spiral structure is that of a Lin-Shu type (Lin 
& Shu 1964), with parameters chosen such that they match 
the RAVE results obtained by Siebert et al. (2012), and it ro¬ 
tates at a constant angular pattern speed that places the sin¬ 
gle co-rotation radius at 12 kpc. However, we note that the 
radial velocity dispersion of the axisymmetric background 
model of F14 is lower than that of the Galaxy, and the radial 
force ratio is high. Therefore, we include another test parti¬ 
cle simulation (Monari et al. in prep.), labelled M15, which 
has a radial force ratio equal to 0.10 and a higher radial 
velocity dispersion. Note that for N-body simulations, the 
spiral arm gravitational potential does not follow a smooth 
cosine curve which is used to describe the spiral arm poten¬ 
tial in the test particle simulations. Therefore, comparison 
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Table 1. Table of simulation parameters. The columns are 1) Model name; 2) Model type; 3) Virial mass; 4) Total disc mass; 5) Radial 
scale length; 6) Spiral arm number; 7) Radial force ratio; and 8) Spiral arm pitch angle. 


Model 

Type 

Mxir (xlQl^Mo) 

Md (xlOiOM©) 

Rd (kpc) 

m 

<? 

7 

K14 

A^-body, co-rotating spirals+bar 

2.5 

5.0 

2.5 

2 

0.16 

10 

F14 

Test particle, DW-like spirals 

0.6 

5.13 

2.0 

2 

0.23 

10 

M15 

Test particle, DW-like spirals 

0.6 

5.13 

2.0 

2 

0.10 

10 

G13F 

A^-body, co-rotating spirals 

1.5 

5.0 

3.5 

4 

0.08 

19 

G14 

A^-body, co-rotating spirals 

1.5 

5.0 

3.5 

6 

0.08 

23 

B15 

Test particle, bar only 

Vc = 220 km s 

fiat rotation curve 

2 

0.015 

- 

Milky Way 

- 

1.0-2.0 

6.0 

2.15 

2-4 

- 

12 


between the q values of A'^-body and test particle simulations 
may not be completely fair. 

The bar model from Bovy et al. (2015), labelled B15, 
consists of a rotating quadrupole (oc cos 2(f>) with a pattern 
speed of 52.25 km kpc“^ and a radial force ratio of 1.5%, 
and makes an angle of 25° with the Sun-Galactic-center 
line. This bar acts as a perturbation to a logarithmic ax- 
isymmetric potential with a circular velocity of 220 kms”*^. 
The bar’s effect on the kinematics of a population of stars 
with a radial velocity dispersion of 31.4kms“^ is computed 
using GALPY^ (Bovy 2015). 


3 THE PECULIAR VELOCITY POWER 
SPECTRUM 

In this section, we outline the procedure for calculating the 
peculiar velocity power spectrum of the simulation data. The 
procedure consists of two steps: 1) to calculate a peculiar ve¬ 
locity field of a given region of a simulated galaxy, and 2) 
to apply a two dimensional Fourier transform to the pecu¬ 
liar velocity field, which is converted to a one dimensional 
power spectrum. Below we give a brief description of the 
implementation of this procedure (see Bovy et al. 2015, for 
more details). 

3.1 The Peculiar Velocity Field 

The total peculiar velocity of a star is defined as the devi¬ 
ation from the mean velocity in the radial, tangential and 
vertical directions, R, 0 and Z, respectively. In a rotation- 
ally supported axisymmetric disc, the mean velocities are 
Vi{ = Vz = 0, and Vg is non-zero and depends on the cir¬ 
cular velocity and the asymmetric drift. We therefore define 
the peculiar velocities of the i-th particle to be 

V_R,i = 

Vg,i = Ve,i-Vg, ( 2 ) 

Vz.i = Vz.i- 

In our coordinate system (with clockwise rotation), positive 
velocities indicate motion away from the Galactic centre in 
the radial direction, motion faster than the mean rotational 
velocity in the tangential direction, and motion in the pos¬ 
itive Z direction. To calculate the LOS components of the 

^ http://github.com/jobovy/galpy . 


intrinsic (or galacto-centric) peculiar velocities defined in 
equation (2), we define an angle, 9, to be the angle between 
the line that connects the star to the galactic centre and the 
line defined by j/ = 0 for positive x, of a cartesian grid whose 
origin is at the galactic centre (see Fig. 1 for a diagram illus¬ 
trating the coordinate system). The peculiar LOS velocity 
can be decomposed into a radial and rotational component 
given by 

'Vio 3 ,r = 'Vr cos bcos{9 - 1), (3a) 

Vios.e = 'Vgcosbsm{9 — 1), (3b) 

where I and b are the Galactic longitude and Galactic lat¬ 
itude coordinates, respectively. Once the radial and rota¬ 
tional LOS peculiar velocities are computed for the star 
particles, the particles are binned onto a grid of spatial res¬ 
olution equal to 0.8 kpc and the mean LOS velocity value is 
taken for each bin. 


3.2 Computing the Power Spectrum 

Once we have calculated the peculiar velocity field, fpq, on 
a grid of N x N points, we calculate the two dimensional 
Fourier transform, Fim, as 

2N 2N 

Firu = Y.Y. (4) 

p =0 ^=0 

where {k^,ky) = {k^^^^l/N, k^°‘^m/N). Here, is de¬ 

fined to be the maximum frequency value that can be sam¬ 
pled (the Nyquist frequency). The grid has been padded at 
the end by N zeros in both directions, in order to remove 
the signal wrap-around pollution caused by the assumption 
of the convolution theorem that the signal is periodic. The 
Fourier-transformed velocity field forms a separate N x N 
grid that stretches from kxy = 0 to fc™’” on both the x 
and y axis, with spacings of Ak = k'^y^ jN. For each value 
Fira = F{kx, ky), the power may be estimated as 

Pikx,ky) = i4nf ■\F{kx,ky)\\ (5) 

The one dimensional power spectrum, P{k), is estimated by 
averaging the power of the two dimensional power spectrum 
in annuli of fc = ^k'^ + ky (see Bovy et al. 2015, for more 
details). 
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Figure 2. Simulation K14: log projected surface density, E (top-left), the magnitude of the intrinsic rotational and radial peculiar velocity 
fields (top-middle and -right panels respectively), the peculiar LOS velocity field (bottom-left), the rotational and radial component of 
the peculiar LOS velocity field (bottom-middle and -right panels respectively). 
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Figure 3. As Fig. 2, but for the spiral analysed in simulation G13F. 
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3.3 The Choice of Solar position 

We apply the methodology described in Sections 3.1 and 3.2 
to snapshots of the simulations described in Section 2. The 
patch of disc in which we choose to apply the Fourier analysis 
extends from —5.5 to 2.5 kpc in X and —3.5 to 4.5 kpc in Y, 
where X and Y represent heliocentric cartesian coordinates 
(as opposed to the galactocentric coordinates, {x, y), referred 
to in the previous section), i.e., (X©, Y©) = (0.0, 0.0). From 
the Solar position, positive X points in the direction of the 
Galactic centre, and positive Y points in the direction of 
rotation (see the dashed box in Fig. 1). The patch provides 
adequate coverage of the spiral arm as viewed from the solar 
position and is similar in size to the disc patch analysed by 
Bovy et al. (2015) and the volume covered by the APOGEE- 
RC data. The bin size of the disc patch is 0.8 kpc, and we 
include all stars within \z\ < 250 pc as in Bovy et al. (2015). 

We set the position of the Sun at the galactocentric 
radius of 8 kpc. The azimuthal angle of the Solar position is 
determined so that there is a spiral arm at a distance of 4 
kpc from the Sun in the direction / = 90 degrees (similar to 
that found by Reid et al. 2014, who fit the spiral arm using 
parallaxes of high mass star forming maser sources). 

We analysed many snapshots in each simulation, and 
found that the power spectrum varies between snapshot for 
the X-body simulations (see Section 4.1). However, the goal 
of this paper is to determine whether a particular model 
can reproduce the features of the observed power spectrum. 
Therefore, for each simulation we focus on the snapshot 
which most closely matches the observational data. Gon- 
trary to the X-body simulations, it is not necessary to ex¬ 
plore many snapshots of the test particle simulations be¬ 
cause the velocity field does not change with time after the 
initial instability caused by the introduction of the rigid spi¬ 
ral potential disappears. Therefore, we show the results at 
a snapshot after the system is well relaxed. 


4 RESULTS 

Fig. 2 shows the peculiar velocity fields and particle distribu¬ 
tion of a snapshot in simulation K14. The plots are centred 
on our chosen Solar position (which is placed at 8 kpc from 
the Galactic centre), and the black box highlights the patch 
of the disc for which we calculate the power spectrum. The 
top-left panel of Fig. 2 shows the log particle surface den¬ 
sity, E. The top-middle and -right panels of Fig. 2 show the 
peculiar rotational and radial velocity respectively, which 
reveals that the spiral structure (as well as the bar) are 
associated with clear features in both velocity fields. For ex¬ 
ample, strong negative radial velocities (towards the galactic 
centre) trace the leading edge of the Perseus-like spiral arm 
located at a distance of about 4 kpc from the solar position 
in the direction I = 90, similar in location to the Perseus 
spiral arm in the Milky Way (Reid et al. 2014), whereas the 
velocity field changes to positive radial velocities (towards 
the galactic anti-centre) on the trailing edge of the spiral 
arm (Grand et al. 2014). 

The bottom-middle and -right panels of Fig. 2 show the 
radial and rotational components of the peculiar LOS ve¬ 
locities. Similar to the intrinsic peculiar velocity fields, the 
peculiar LOS velocity fields show significant fluctuations, al¬ 
though the fluctuations are significantly reduced in regions 


where the LOS velocity is not sensitive to the rotational and 
radial velocity, i.e., along the X-axis for Y — 0 and an arc 
region whose galactocentric radius is similar to that of the 
Sun, respectively. The bottom-left panel of Fig. 2 shows the 
total LOS peculiar velocity that results from the summation 
of the rotational and radial LOS components. Figs. 3, 4, 5 
and 6 show the same velocity fields as Fig. 2 for the spirals 
examined in G13F, F14, M14 and G14, respectively. All sim¬ 
ulations show velocity fields that trace the spiral structure. 
Similar to K14, the radial and rotational components of the 
intrinsic peculiar velocity fields of G13F (Fig. 3) and G14 
(Fig. 6) exhibit radial velocities that generally point out¬ 
ward (inward) on the trailing (leading) side of spiral arms, 
and rotational velocities which are slower (faster) than the 
mean on the trailing (leading) side of spiral arms (Grand 
et al. 2014). These trends of intrinsic peculiar velocities is 
evident in the LOS velocities as well. The main difference 
that can be seen by eye between Figs. 2 and 3 is the lack of 
velocity structure in the central region of the galaxy (where 
4 < X < 8 and —4 < X < 4) owing to the presence (lack 
of) a bar in K14 (G13F). In comparison to the X-body sim¬ 
ulations, the test particle simulations F14 (Fig. 4) and M15 
(Fig. 5) show very clear spiral patterns in the surface density 
and velocity fields, which arises from the analytically defined 
spiral density perturbation. In these models, the mean ra¬ 
dial velocity inside the spiral arms generally points inward 
(negative velocity), whereas it points outward (positive ve¬ 
locities) in the inter arm regions (Siebert et al. 2012). The 
main difference between the peculiar velocity fields of F14 
and M15 is that the magnitude of the fluctuations is larger 
for F14, in particular in the rotational component, which 
is a result of the larger spiral amplitude of F14 relative to 
M15. 

Fig. 7 shows the azimuthally averaged one-dimensional 
power spectra of the velocity fields for K14 shown in Fig. 
2, in addition to the power spectra of the velocity fields in 
a similar region of the other simulations. The top-left panel 
plots power spectra of the simulations alongside the best 
fitting bar model of Bovy et al. (2015), B15, and spectra 
calculated from the observed peculiar velocities from the 
APOGEE and RAVE surveys (Bovy et al. 2015). The power 
spectrum of the fiducial simulation K14 exhibits a turnover 
at peak power in the range 0.3 - 0.4 kpc”*^, and is a good 
match to the data and bar model of B15, although the power 
at scales of fc ~ 0.3 kpc“^ appears to be larger than that 
of the data by about 2 km s“^. We have confirmed that 
the power spectrum does not change significantly for the 
case of a more limited volume like the region covered by the 
APOGEE-RG data. 

Interestingly, the power spectrum of the peculiar LOS 
velocity field in the spiral-only N-body simulation, G13F, 
exhibits a similar peak shape to the observed power spec¬ 
trum, although the amplitude on all scales is smaller than 
that of K14. This reflects the strength of non-axisymmetric 
structures, which is larger for K14 than for G13F (see the q 
column of Table. 1). In contrast to the X-body simulations, 
the power spectra of the density wave models of F14 and 
M15 show increasing power on progressively larger scales. 
This indicates that the density wave-like spiral models do 
not fit the data well, which is consistent with other density 
wave - test particle models shown in Bovy et al. (2015). 

The above results show that in addition to the bar only 
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Figure 4. As Fig. 2, but for the spiral analysed in simulation F14. 
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Figure 5. As Fig. 2, but for the spiral analysed in simulation M15. 
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Figure 6. As Fig. 2, but for the spiral analysed in simulation G14. 


model, B15, the barred-spiral A'^-body simulation, K14, is 
able to reproduce the observed peculiar LOS power spec¬ 
trum very well, in contrast to the cosmological Af-body sim¬ 
ulation (Guedes et al. 2011) investigated by Bovy et al. 
(2015). This suggests that the observed peculiar motions can 
be explained by either a combination of a bar and transient, 
co-rotating spiral arms (Grand et al. 2012b; Kawata et al. 
2014) or the bar component only as found by (Bovy et al. 
2015, B15). In contrast, the observed power spectrum cannot 
be explained by the stellar motion induced by the density 
wave-like spiral arms in test particle simulations F14 and 
M15. On the other hand, the power spectrum of simulation 
G13F (and to an extent G14) shows evidence of a turnover 
in power on scales of A: ~ 0.4 km s~^, which is the charac¬ 
teristic feature of the observed power spectrum. Therefore, 
it is possible that the main features of the observed power 
spectrum can be explained purely by the systematic stel¬ 
lar motion induced by co-rotating spiral arms seen in sim¬ 
ulations (e.g.. Grand et al. 2012a,b, 2014). We note that 
this study is limited to our current set of simulations, and 
we have not extensively searched for a Af-body simulation 
model with transient spiral arms only that can better repro¬ 
duce the observed power spectrum. By the same token, we 
have analysed only two test particle simulations with den¬ 
sity wave-like spiral arms, and therefore cannot state out¬ 
right that no density wave configuration can explain the 
observed power spectrum. Our study highlights that the ob¬ 
served power spectrum of the APOGEE-RC data can be 
explained not only by the bar induced velocity fields, but 
also by the velocity field induced by the spiral arms. It is 
therefore interesting to explore many more different models 
of spiral arms and compare them with the observations. 


The middle (bottom) panels of Fig. 7 show the power 
spectra of the radial and tangential components of the LOS 
(intrinsic) peculiar velocity fields. The power spectra of 
model F14 and M15 show much less power in the intrin¬ 
sic velocity fields than the simulation K14 and G13F, which 
results in less power in the LOS rotational and radial ve¬ 
locities, the latter of which seems to drive the shape of the 
power spectrum in the top-left panel. 


It is interesting to note that the power in K14 is dom¬ 
inated by the radial LOS velocity, despite the relatively 
strong intrinsic rotation peculiar velocity. This severe reduc¬ 
tion in power between intrinsic and LOS rotational velocity 
may be explained by the ‘flip’ in LOS velocity direction that 
occurs either side of T = 0 in the bottom-middle panel of 
Fig. 2. For example, the positive intrinsic rotation velocity 
region at {X, Y) = (2, 4) in the top-middle panel of Fig. 
2 corresponds to a receding (positive) LOS velocity in the 
bottom-middle panel of Fig. 2. However, the strong nega¬ 
tive intrinsic rotation velocity region at {X, Y) = (—2, —2) 
in the top-middle panel of Fig. 2 indicates a slow rota¬ 
tional peculiar velocity, which leads to a receding (positive) 
LOS velocity as viewed from the solar position. It appears 
that this ‘flip’ in velocity direction caused by the conver¬ 
sion from intrinsic to LOS velocity can affect the resultant 
power spectrum, which therefore is dependent on the se¬ 
lected area of analysis and does not reflect the true effects 
of non-axisymmetric structure on the intrinsic radial and ro¬ 
tational velocity fields. It is therefore important to observe 
the 3D velocity held in a larger area, which will soon be 
possible for the observational data provided by Gaia. 
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Figure 7. The azimuthally averaged one-dimensional power spectrum of: the peculiar LOS velocity field (top-left), the rotational and 
radial components of the peculiar LOS velocity field (middle-left and -right) and the intrinsic rotational and radial peculiar velocities 
(bottom-left and -right). Each solid curve denotes the power spectrum of the respective field from a simulation model. The observed 
power spectrum measured by Bovy et al. (2015) from APOGEE-RC and RAVE data is shown in black points with associated error bars. 
The shaded area marks the 95% noise region of the data. 


4.1 Other Power Spectrum dependencies 

4.1.1 Spiral arm number 

Fig. 7 contains the power spectra of simulation G14, which 
is a flocculent m = 6 spiral galaxy simulation suitable for 
comparison with G13F and K14 in order to study the ef¬ 


fect of spiral arm number on the power spectrum. The most 
pronounced difference is visible in the bottom-right panel 
of Fig. 7, which shows that the power in intrinsic peculiar 
radial velocity peaks at smaller scales, k = 0.5 kpc“^, than 
G13F and K14, which peak k — 0.3 - 0.4 kpc“^. This in¬ 
dicates that flocculent spiral structure produces more power 
on small scales owing to the smaller inter arm separation. 
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Figure 8. The same as Fig. 7 but comparing the power spectra of 
the fiducial spiral arm of K14 from two different solar positions: 
a position nearer to (blue) and a farther from (red) the spiral in 
comparison to the fiducial K14 position shown in Fig. 7. 


which limits the spatial extent of the peculiar velocity fluc¬ 
tuations. However, the power spectrum of the total LOS 
peculiar velocity is not distinguishable in shape from simu¬ 
lations K14 and G13F. 


4-1.8 Proximity of Spiral arm 

Fig. 8 shows how spiral arm proximity with respect to the so¬ 
lar position affects the power spectrum in simulation K14. 
We consider two conhgurations in addition to the fiducial 
power spectrum: one in which the disc is rotated 20 de¬ 
grees clockwise from the hducial setup, and one in which 
is rotated 20 degrees anti-clockwise from the fiducial setup. 
In this way we obtain three solar positions with respect to 
the spiral arm: a far, an intermediate (fiducial) and a near 
position. It is clear from Fig. 8 that neither the near nor 
far conhguration reproduces the observed power spectrum; 
the far configuration, though it reveals similar intrinsic pe¬ 
culiar velocity power spectra, does not reach the height of 
the observed peak of the LOS peculiar velocity spectrum, 
whereas the near configuration exhibits a contrasting power 
spectrum in all peculiar velocity fields. Of particular note is 
the double peak structure in radial and total LOS peculiar 
velocity, which may arise from the faint stellar filamentary 
structure on the leading side of the main ‘Perseus’-like spi- 



Figure 9. The same as Fig. 7 but comparing the power spectrum 
of spiral arm G14b. The blue corresponds to the spiral arm at the 
earlier stages of formation when the pitch angle is higher, whereas 
the red line corresponds to the same spiral arm at a later stage 
of formation when the pitch angle is lower. 


ral. Such a structure in close proximity of the main spiral 
may produce similar effects to the flocculent spiral structure 
described above. 


4-1.3 Evolutionary stage of the Spiral arm 

As highlighted by Baba et al. (2013) and Grand et al. (2013), 
spiral arms in A’-body simulations continuously wind up 
from the moment that they begin to form until they are 
disrupted, and therefore each spiral arm sweeps though a 
range of pitch angles during its lifetime. In this section, we 
explore the dependency of the power spectrum on different 
spiral evolutionary stages (and therefore, pitch angle), which 
has been previously reported to affect the velocity distribu¬ 
tion of the solar neighbourhood (De Simone et al. 2004). We 
focus on a spiral arm from simulation G14 whose winding 
evolution is studied in detail by Grand et al. (2014). Note 
that this spiral arm is a different spiral arm from the spiral 
studied above (labelled G14), therefore we denote this spiral 
arm G14b to distinguish it from the original. Spiral G14b is 
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in a less flocculent region of the disc^, which facilitates the 
stndy of the effect of pitch angle. 

Fig. 9 shows the power spectra of the same size velocity 
field patch at two different evolutionary stages of the same 
spiral arm: an early stage during which the spiral has a rela¬ 
tively high pitch angle, and a later stage when the spiral arm 
is relatively tightly wound. Both the early and late stage spi¬ 
ral arms produce very similar LOS velocity power spectra. 
However the intrinsic radial and rotational velocity fields 
reveal that the spiral of larger pitch angle (G14b-early) pro¬ 
duces more power in the radial peculiar velocity field than 
the spiral of lower pitch angle at the later spiral evolntion- 
ary stage (G14b-late). It seems plansible that there exist 
spiral-driven time-dependent peculiar velocity fields, which 
evolve with the winding spiral. Indeed, Fig. 10 demonstrates 
a scenario in which the open spiral arm (G14b-early) drives 
strong systematic radial motions over a large radial range, 
which is reflected in the bottom-right panel of Fig. 9. As the 
spiral arm winds to lower pitch angles (G14b-late), the ve¬ 
locity vectors appear to re-orient themselves to a more tan¬ 
gential direction, which has the effect of reducing the power 
in the radial peculiar velocity spectrum while increasing the 
power in the rotational peculiar velocity spectrum (though 
the difference shown in the bottom-left panel of Fig. 9 is 
slight). The leading cause of this difference may be related 
to either the different stages of spiral arm formation and dis¬ 
ruption (see, for example. Baba et al. 2013), or the different 
pitch angle of the spiral arm at these instances. A further 
detailed study is required to determine the cause. 

Although the above analysis indicates that the peculiar 
velocity field is dependent on spiral arm number, solar po¬ 
sition and spiral arm pitch angle, we caution that further 
work covering a much larger suite of simulations is required 
to fully work out the dependencies and degeneracies of the 
peculiar velocity power spectrum. On the other hand, this 
means that the stage of spiral arm evolution may be able 
to be constrained by analysis of the peculiar velocity fields 
presented in this paper. 

4.2 External galaxies 

The power spectrum analysis presented so far is a useful 
tool for quantifying the perturbed velocity field in the Milky 
Way, in which we can observe only a local patch of the disc 
at present (see also Bovy et al. 2015). For external galaxies, 
however, we have available the whole disc region on which to 
apply the power spectrum analysis. Therefore different re¬ 
gions of the disc, such as the centre, can be examined as well 
as the whole disc. The power spectrum analysis of a region 
that extends over the whole disc is more robust than that of 
the local disc patch analysed in the previous sections, which 
can be sensitive to the choice of patch (see Fig. 8). In Fig. 
11, we show the power spectra of the rotational and radial 
peculiar velocity fields for different regions of the disc. In 
the top panel of Fig. 11, we present the power spectra of the 
central 8x8 kpc region of the disc, i.e, extending from 4.0 
to 12.0 kpc in X and from —4.0 to 4.0 kpc in Y, for models 

^ Although simulation G14 (Grand et al. 2014) is reported to 
have m = 6 spiral arms, this number is not always constant with 
time nor is the spiral structure perfectly symmetric. 



e 


Figure 10. Close up density map in polar coordinates of the spi¬ 
ral arm in simulation snapshots G14e (upper) and G141 (lower). 
The top panel shows the forming spiral arm in front of the solar 
position at {R,6) ~ (8,150). The bottom panel shows the same 
spiral arm at a later time when it is fully formed. Over-plotted is 
the mean peculiar velocity field, which is shown to change orien¬ 
tation as the spiral arm evolves. 


K14 and G13F (we do not show the power spectra of the 
density wave models, F14 and M15, because these models 
follow particle motion at i? > 4 kpc only). For K14 we anal¬ 
yse also the power spectrum of the gas velocity field. There 
are two main characteristics of the peculiar velocity power 
spectrum of the central region: the first is the large power in 
the fluctuations of both rotational and radial peculiar veloc¬ 
ity fields, and the second is that the radial peculiar velocity 
fluctuates on the largest scales of ~ 0.2 kpc”*^, whereas the 
rotational peculiar velocity field peaks around the smaller 
scale regime of fc ~ 0.5 kpc“^. The gas shows even more 
pronounced velocity fluctuations. The power spectra for the 
bar only model, B15, is shown to have a similar shape to 
those of K14, which is expected owing to the presence of 
a bar in the central regions of both models (although we 
caution that the bar in B15 is modelled only as a rotating 
quadrupole potential oc cos 20 in the azimuthal direction). 
In contrast, model G13F shows little power, given the lack 
of non-axisymmetric structure in the centre. 

The middle panel of Fig. 11 shows the power spectra of 
the peculiar velocity fields for a region that we define as the 
‘whole disc’, i.e., extending from -4.0 to 20.0 in X and from 
-12.0 to 12.0 in Y. The relatively large size of the ‘whole disc’ 
patch makes it possible to probe scales up to A: ~ 0.05 kpc“^. 
For our fiducial spiral model, K14, the power decreases on 
all scales in comparison to that of the central region, and 
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Figure 11. The power spectrum of the intrinsic peculiar rota¬ 
tional and radial velocity fields of the central region of the disc 
(top), the whole disc (middle) and the whole disc with the central 
region cut out (bottom). Note the change in the x-axis and y-axis 
ranges. 


the shape of the radial peculiar velocity power spectrum 
shows a single prominent peak at around scales of fe = 0.2 
kpc“^, with decreasing power on scales k < 0.1 kpc“^. For 
the G13F spiral, the power of peculiar velocity fields is again 
lower than that of K14, particularly the radial velocity field, 
though the spectra exhibit a broad peak shape which is a 
consequence of the spiral structure. The power spectra of 
both fields are peaked at scales of around k ~ 0.2 kpc“^, 
although the maximum power is only 5 km s“^, much less 
than the 18-20 km s“^ peak of K14. For B15, the distinct 
peak observed in the top-left panel is no longer present, while 
the power spectrum of the radial peculiar velocity field is 
reduced and distinguishable from that of K14. 

The bottom panel of Fig. 11 shows the power spectra 


for the whole disc region (as described above), in which the 
peculiar velocities in the central region (the same central 
region analysed in the top panels) are set to zero. We have 
re-introduced the density wave spiral models as these now 
contain the required data. The shape of the power spec¬ 
tra for all models does not appear to change from that of 
the whole disc upon the exclusion of the central region, and 
shows reduced power in most cases. The rotational peculiar 
velocity power spectra reveal that all simulations, except 
the gas component of K14, show similar features of increas¬ 
ing power with larger scales, and are therefore difficult to 
distinguish. 

The radial peculiar velocity power spectrum yields 
striking differences between the modefs presented; the 
barred-spiral model of K14 shows a highly peaked power 
spectrum that peaks on scales oi k — 0.2 kpc“^, whereas 
model G13F shows a broader peak that reaches a maximum 
on similar scales. In contrast, models F14 and M15 exhibit 
no such peak in their power spectra, with very little power 
on all scales. This analysis indicates that the combination 
of both radial and rotational peculiar velocity fields in var¬ 
ious regions of the discs of external galaxies is required to 
distinguish between different spiral arm and bar models. 


5 CONCLUSIONS 

We have presented a quantitative analysis of the scale and 
power of the fluctuations in peculiar velocity induced by 
non-axisymmetric structure, which we apply to a suite of 
numerical simulations. The suite comprises N-body simula¬ 
tions that naturally produce transient, winding spiral arms 
of varying types, and test particle simulations that treat spi¬ 
ral structure as an analytically prescribed density wave. We 
compare the results to the latest Milky Way data, and come 
to the following conclusions: 

• From the fiducial Milky Way-type N-body simulation 
K14, we select a spiral arm that is similar to the Perseus 
arm of the Milky Way, and determine the solar position such 
that the spiral arm is ~ 4 kpc from the Sun in the direction 
I — 90, similar to what is found by Reid et al. (2014). We 
show that it is possible for the combination of a bar and 
transient, co-rotating spiral arms to reproduce the observed 
power spectrum of the APOGEE-RC and RAVE data. How¬ 
ever, the test particle simulations with density wave-like spi¬ 
rals do not fit the observed power spectrum well, which is 
consistent with previous work (Bovy et al. 2015). However, 
it remains to be seen whether or not the effects of density 
wave spiral arms combined with a bar can reproduce the 
data. Therefore, analysis of test particle simulations that 
include a bar and spiral wave potential are a natural next 
step to this work, and will be presented in a forthcoming 
paper (Monari et al. in prep.). 

• It is possible that transient, winding spiral structure 
alone can reproduce the defining peak feature of the ob¬ 
served power spectrum. This emphasises the differences be¬ 
tween the effects of transient co-rotating spiral arms and 
density waves, and suggests that the power spectrum of K14 
is not purely caused by the bar. 

• The peculiar velocity power spectrum is sensitive to the 
proximity of the spiral arm with respect to the solar position, 
the number of spiral arms (or the degree of flocculence) and 
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the pitch angle of the spiral arm. However, further study is 
needed to separate out all of these dependencies fully. 

• Analysis of the radial and rotational peculiar velocity 
fields in the central regions of discs shows that barred models 
can easily be separated from spiral only models. 

• The power spectra of whole disc regions are a more 
robust way of characterising the velocity fluctuations, and 
appear usable to distinguish between the simulations pre¬ 
sented in this paper. Furthermore, gas motions are shown 
to follow those of the stars and exhibit very similar power 
spectra, though with more power in general. 

In this paper, we have calculated the power spectra of 
Fourier-transformed peculiar velocity fields primarily as a 
means to quantify and compare with the observational and 
model data of Bovy et al. (2015). However, we have shown 
that information of both the intrinsic radial and tangen¬ 
tial peculiar velocity field is required in order to distinguish 
between the different mechanisms to induce the peculiar ve¬ 
locity field in the Galactic disc. Therefore, it is necessary for 
Galactic surveys to go beyond the LOS velocity informa¬ 
tion. The full phase space information that surveys such as 
Gaia will yield is expected to provide excellent constraints 
on dynamical models of the Galaxy and non-axisymmetric 
structure thereof. 

In addition to the Milky Way, we have also shown that 
we can make predictions for external galaxies, for some of 
which peculiar velocity information should be readily avail¬ 
able from integral field spectroscopic surveys such as MUSE. 
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